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We consider the T- matrix near the unitarity limit and the energy dependence of the saturation 
momentum. We discuss the solution to the Kovchegov equation, or equivalently, to the BFKL 
evolution in the presence of a single saturation boundary. We include some of the correlations 
missed in the Kovchegov equation by solving the BFKL equation in the presence of two 
boundaries. The T-matrix now turns out to be frame-independent, which was not the case 
for the solution in the case of a single boundary, and it doesn't show the scaling behavior 
of the solution to the Kovchegov equation. We find for the saturation momentum an energy 
dependence which differs from the one following from the Kovchegov equation. 



1 Introduction 

The growth of cross sections with increasing energy, or parton densities with decreasing Bjorken- 
x, at a fixed hard scale is given by the BFKL^ evolution when parton occupation numbers are 
not too large or the T-matrix not to close to the unitarity limit. At or near the unitarity 
limit nonlinear parton evolution becomes important which is not described by the linear BFKL 
equation. An extension of the BFKL equation that includes also nonlinear evolution is the 
Balitsky equation^ or the Jalilian-Marian, Iancu, McLerran, Leonidov and Kovner (JIMWLK)^ 
equation. The Balitsky and JIMWLK equations are coupled equations involving higher and 
higher correlations and as such are very difficult to deal with analytically. A better understanding 
of the Balitsky and JIMWLK evolution may emerge from numerical calculations such as^l. 

Kovchegov^ has suggested a simpler equation than the Balitsky or JIMWLK equations to 
deal with scattering at or near the unitarity limit. The Kovchegov equation can be viewed 
as a mean field version of the Balitsky or JIMWLK equation in which higher correlations are 
neglected. While incomplete, as is any mean field like approximation, the Kovchegov equation 
is likely the best equation one can write down in terms of functions which has built in correct 
unitarity limits for high energy scattering. Many interesting limits of the Kovchegov equation 
have been understood by analytical methods. In this work we focus on the transition region from 



weak to the saturation regime in which Munier and Peschanksi - ' have determined the form of the 
T-matrix and the rapidity dependence of the saturation momentum by solving the Kovchegov 
equation analytically. The same results have been obtained even before by the authors of Ref.^ 
by solving the BFKL equation in the presence of a single saturation boundary which effectively 
approximates the unitarity limit guaranteed in the Kovchegov equation. In the next Section we 
derive the Kovchegov equation and study its limitations due to the missed higher correlations. 

We attempt to go beyond the mean field like approximat ion in the Kovchegov equation by 
solving the BFKL equation in the presence of two boundaries'^. We find in the vicinity of the 
saturation regime an expression for the T-matrix which is frame-independent and doesn't show a 
scaling behavior in contrast to the solution to the Kovchegov equation. The energy dependence 
of the saturation momentum is also different from the one coming from the Kovchegov equation. 



2 The Kovchegov equation 

Consider the high-energy scattering of a color dipole on a target (another dipole, hadron, or 
nucleus) at relative rapidity Y in a frame where the dipole is an elementary quark-antiquark 
pair and the target a highly evolved system. Now we wish to know how the elastic scattering 
amplitude S(x j0 ,x_i,Y) changes when the rapidity Y is increased by a small amount dY (x 
and x_i are the transverse coordinates of the quark and antiquark of the dipole) . The change of 
S(x ,x_i,Y) with rapidity is determined by the Balitsky or JIMWLK equation 



^S(x 01 ,Y) = ^J d 2 x 2 ^- [s^(x 02 ,x 12 ,Y) - S(x 01 ,Y)} . (1) 



This equation can be interpreted as follows: When the increase dY is viewed as increasing the 
rapidity of the dipole then the probability for the dipole to emit a gluon (transverse coordinate 
x 2 ) increases. In the large N c limit this quark-antiquark-gluon state can be viewed as a system 
of two dipoles. The scattering of this two dipole state on the target is given by S^ 2 \x J02 , x 12 , Y) 
whi le th e subtracted S(x i,Y) is the virtual contribution necessary to normalize the wavefunc- 
tionEl] The later gives the scattering of a single dipole on the target because the gluon is not 
in the wavefunction of the dipole at the time of the scattering. The weig ht in Eq. □ gives the 
probability for the emission of two dipoles of and x_i 2 from the initial dipole of size Xqi . 

It is difficult to use the Balitsky- JIMWLK equation because of the unknown (x 02 ,x 12 , Y). 
The assumption that the scattering of the two dipole state on the target factorizes 

S^(x 02 ,x 12 ,Y) = S(x 02 ,Y)S(x 12 ,Y) , (2) 

which is a sort of a mean field approximation for the gluonic fields in the target, leads to the 
Kovchegov equation^ 



^S(x 01 ,Y) = ^J d2 ^^T [S(x 02 ,Y)S(x 12 ,Y) - S(x 01 ,Y)} . (3) 



One exact result of the Kovchegov equation for the 5-matrix deep in the saturation region 
is the Levin- Tuchin formula^ 

S(p,b,Y) ~ e ~ c(p - ps)2 (4) 

im 

with the constant c = — Cp(l — Ao)/(iV c 2x(Ao)). Recently, the authors of Ref. - have claimed 
that S deep in the saturation regime has the form given by Eq. 0] but with a constant at least 
a factor of 2 smaller than the c which follows from the Kovchegov equation. The cause for this 
discrepancy is the lack of fluctuations in the Kovchegov equation. 

Now let us focus on the region close to the saturation regime. Consider the scattering of 
a dipole of size x on a dipole of size x' at relative rapidity Y and impact parameter b. Near 



the saturation boundary the Kovchegov equation or the BFKL equation in the presence of an 
absorbtive saturation boundary give for the T-matrix in laboratory frame 



T{x,x',Y;b) =T (b,x) Q 2 S (Y) x 
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and for the rapidity dependence of the saturation momentum 
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where Ao = 0.372. Near the saturation boundary, i.e., for the size of the probe x close to 
1/Q 2 (Y), the scattering amplitude in © shows a scaling behavior since it depends on Q 2 S {Y) x' 2 ; 
lines with constant Q 2 (Y) x' 2 are lines of constant scattering amplitude. 

Recently^^ we have shown that fluctuations are important in evolution also in the region 
near the saturation boundary ("scaling region"). We have found that the completeness relation 



is not satisfied for the dipole number density n which is obtained by using the BFKL evolution 
in the presence of a single saturation boundary, as in the case of the T-matrix in Eq. This 
mismatch is equivalent to the statement that the result for the T-matrix in different frames 
is different. The reason for this mismatch is the mean field like treatment of fluctuations in 
the Kovchegov equation: On the left hand side of Eq. the evolution from the initial point 
(x, y = 0) to the final point (V, y = Y) occurs in one step while on the right hand side of Eq. [7| 
the evolution can be viewed as proceeding in two steps, form (x,y = 0) to (r2,Y/2) to (x',Y), 
and at the intermediate stage of evolution the fluctuations are not properly taken into account. 

In terms of BFKL evolution in the presence of a saturation boundary the reason why the 
completeness relation is not satisfied goes as follows: on the lhs of Eq. the entire evolution 
occurs in the presence of a fixed boundary Q s determined by x and y (Fig. 0t) while on the 
rhs of Eq. the second part of evolution occurs in the presence of a different boundary Q s 
determined by r 2 and y (Fig. 0)). Since the dipole number density is obtained in the mean field 
approximation, one expects that the completeness relation should be fulfilled in an approximative 
way. Indeed, if one requires for the evolution on both sides of Eq. 0to happen in the presence 
the same saturation boundary (Fig. ^i), i.e., if one extends the diffusion region for the second 
part of the evolution on the rhs of Eq. from Q s to Q s , then the lhs and rhs of Eq. give the 
same result. However, because of the extended diffusion region, there are now paths of ev oluti on 
from (r 2 ,Y/2) to (x',Y) which manifestly violate unitarity (T ^> 1). See for more details^^. 

3 BFKL evolution in the presence of two saturation boundaries 

Now we introduce a second absorptive barrier at P2(y), as shown in Fig. such that, all 
unitarity violating evolution between the initial point (x, 0) and the final point (x', Y) is 
eliminated when the evolution is viewed as proceeding in two steps The boundary p\(y) 
corresponds to saturation in the wavefunction of the evolved dipole x scattering on an elementary 
dipole x' (forward evolution) and P2(y) corresponds to saturation in the wavefunction of the 
evolved dipole x' scattering on the elementary dipole x (backward evolution). The introduction 
of the second boundary p2{y) and the symmetry it brings with has another benefit: The T-matrix 
is now frame-independent which was not the case for the single boundary case. 



n(x,x',Y) = — i n(x,r 2 ,Y/2) n(r 2 ,x',Y/2) , 
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Figure 1: BFKL evolution in the presence of shaded saturation regions in the Y — ln(x 2 /z 2 ) plane with z generic: 
(a) evolution in the presence of a single saturation boundary, (b) the second evolution step from Y/2 — > Y has its 
own saturation boundary, (c) evolution in the presence of two boundaries. 



Nea r the saturation regime the BFKL evolution in the presence of two absorbtive bound- 
aries^^ gives a T-matrix which depends on ]n(Q s (Y) x')/[aY/(Ap) 3 } with Ap ~ i-\ hi(l/a). 
Thus, the scaling behavior of the solution to the Kovchegov equation is now lost. The rapidity 
dependence of the saturation momentum for aY 3> (p2 — pi) 2 / (irN c x" (A^)) becomes 

Q 2 s (x, Y; b) = Qq(x, b) \ exp 
with Arf = Aq + 7r 2 /[2(A/9) 2 (l — Aq)] and is different as compared to the result in Eq. HO 
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